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ABSTRACT 

Aims. Dust properties appear to vary according to the environment in which the dust evolves. Previous observational indications of 
these variations in the far-infrared (FIR) and submillimeter (submm) spectral range are scarce and limited to specific regions of the 
sky. To determine whether these results can be generalised to larger scales, we study the evolution in dust emissivities from the FIR 
to millimeter (mm) wavelengths, in the atomic and molecular interstellar medium (ISM), along the Galactic plane towards the outer 
Galaxy. 

Methods. We correlate the dust FIR to mm emission with the HI and CO emission, which are taken to trace the atomic and molecular 
phases, respectively. The study is carried out using the DIRBE data from 100 to 240 yum, the Archeops data from 550 fim to 2.1 mm, 
and the WMAP data at 3.2 mm (W band), in regions with Galactic latitude |b| < 30°, over the Galactic longitude range (75° < 1 < 198°) 
observed with Archeops. 

Results. In all regions studied, the emissivity spectra in both the atomic and molecular phases are steeper in the FIR (6 = 2.4) than 
in the submm and mm (fi = 1.5). We find significant variations in the spectral shape of the dust emissivity as a function of the dust 
temperature in the molecular phase. Regions of similar dust temperature in the molecular and atomic gas exhibit similar emissivity 
spectra. Regions where the dust is significantly colder in the molecular phase show a significant increase in emissivity for the range 
100 - 550 f^m. This result supports the hypothesis of grain coagulation in these regions, confirming results obtained over small 
fractions of the sky in previous studies and allowing us to expand these results to the cold molecular environments in general of the 
outer MW. We note that it is the first time that these effects have been demonstrated by direct measurement of the emissivity, while 
previous studies were based only on thermal arguments. 

Key words. Dust, extinction - Infrared: ISM - Submillimeter 



1. Introduction 



Measuring the dust emissivity is important to inferring the nature 
of dust from its emission and also determining the dust heating 
from its observed temperature. In addition, variations in the dust 
emissivity may seriously affect the mass estimates inferred from 
FIR to mm observations. Dust emissivity and its possible vari- 
ations with wavelength or temperature are also critical for sep- 
arating astrophysical foreground emission from the cosmic mi- 
crowave background (CMB). However, deriving dust emissivity 
is generally a difficult task because it relies on both estimating 
the dust temperature and the gas column den sity associated with 
the emitting dust. Big dust Grains (BG, see iDesert et al.L fl990. 
for a description of each dust component) are in thermal equi- 
librium with the interstellar radiation field (ISRF). Their emis- 
sion is close to that of a gray-body wi th an equilibrium tem- 
perature near 1 7.5 K in the diffuse ISM dBoulanger et al.Lll996t 
lLagache et all Il998l) . with a maximum in the far-infrared. The 
emission spectrum, assuming a fixed dust abundance and a sin- 
gle grain size, follows 

I V (A) = e v (A)B v (A, T d )N H , (1) 

where \ V { A) is the specific intensity or brightness (energy flux per 
unit frequency and solid angle), B v is the Planck function, Tj is 
the equilibrium temperature of grains, Nh is the total column 
density (number of hydrogen atoms per unit area), and e v (A) is 
the wavelength-dependent dust grain emissivity. 



Although, in principle, t he emissivity may depend on wave- 
length and temperature (e.g., Me ny et aUl2007l) . it has been cus- 
tomary to assume no temperature dependence and a power law 
distribution with frequency, 



(2) 



€ V (A) = X d €o \ — 



where is the dust to gas mass ratio, e () is the emissivity at 
wavelength Aq, and B is the emissivity spectral index. The spec- 
tral index is usually assumed to equal 2, which is correct for 
crystalline grains. Howeve r, dust grains are kno wn to be amor- 
phous in the diffuse ISM dKemper et al.L 12004 and this could 
have an impact on the emissivity. Grains are probably affected 
by different physical processes depending on the environment. 
Their properties are likely to vary according to the density of the 
medium. Different studies have highlighted changes in the grain 
properties in the den se and cold medium , such as the depletion o f 
chemical elements (Ca selli et al.L 119991 : iBacmann et all 12002). 
that imply that condensation of gaseous species is occuring on 
the grain surface. 

A decrease in the IRAS I60A100 rati o in several isolated 
molecular clouds has also been noted (e.g jLaureijs et all [l 99 1 ; 
lAbergel et all 11994 1 19951) . whereas this ratio seems to be con- 
stant in diffuse regions. This could be interpreted as a decrease 
of very small grain (VSG) abundance with respect to that of the 
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BGs in dense molecular clouds. Bern ard et al.l (1999) analysed 
the FIR dust emission from a molecular cirrus in Polaris, using 
the PRONAOS balloon-borne experiment data and found a larger 
decrease in the BG equilibrium temperature than predicted from 
the decrease in the ISRF intensity in the cloud. The same region 
showed an obvious deficit in 60 /im emission. They attributed 
the change in the particular dust properties of dense environ- 
ments; 

IStepnik et al.l d2003l) attempted to explain the origin of the 
decrease in VSG abundance and the drop in the BG temperature 
in dense environments. They measured a BG temperature of 16.8 
K outside a filament in Taurus, and of 12 K at the centre of the 
filament, using the PRONAOS and the IRAS data. They showed 
that a model that assumed the standard properties of dust, with- 
out any spatial variations in the grain properties, was unable to 
reproduce the observations. The attenuation of the radiation field 
in the cloud was also unable to explain the shape of the observed 
emission profiles. They proposed that the aggregation of large 
dust grains and most of the VSGs in fractal assemblies would 
explain both the observed 60 /im emission deficit and the unusu- 
ally cold BG temperature. They showed that aggregation would 
lead to an increase in the BG emissivities by a factor of 3 - 4, 
and that between 80 and 100% of the VSGs should participate 
in the aggregation process. The best-fit model showed that the 
VSG and BG aggregation occurs at an extinction Av > 2.1 mag 
and at densities higher than ne=3 H/cm 3 . However, the emissiv- 
ity increase in dense environments remains indirect evidence of 
change in the emissivity, since Tj is lower than expected for the 
corresponding ISRF. This effect has never been illustrated by a 
direct determination of the dust emissivity, mainly because both 
the dust t emperature and Nh are difficult to measure. 

Camb risy et al.ld2001al) analysed the DIRBE emission of the 
Polaris complex over 20 deg2, to determine the dust properties. 
They decomposed the emission into a warm and a cold compo- 
nent constrained using the I60A100 IRAS intensity ratio. They de- 
termined the extinction map o f this complex in the visible, by us- 
ing star counts in the V band ( Cambres vl 1 1 9991) . They calculated 
the submillimeter emissivity of each component, normalised to 
the visible extinction, and found that the emissivity is 4 times 
higher in the cold region of the complex than in the warm re- 
gion. This result matches the emi ssivity increase deriv ed to ex- 
plain the low dust temperature by Step nik et al.1 (120031) . and was 
based on the determination of the submm to visible dust opacity 
ratio, which traces the optical properties of the grain material, 
independent of the gas column density. 

Two hypotheses have been prop osed to explain the observa- 
tions: the formation of ice mantles ()Laureiis et all 1 1 99 ll T l 996 ) 
and/o r grain aggregation dDraine & Anderson! 119851; iTielensl 
1989). The 60 /im emission deficit in dense regions can be at- 
tributed to a large majority of VSGs participating in the aggre- 
gation process, and the VSGs are almost totally accreted onto 
the ag gregates, in the inner regions of the cloud (Stepnik et al., 
2003). The VSGs are then connected to the aggregate, which 
acts as a thermostat. Their temperature then stops fluctuating 
and their contribution to the mid-IR greatly decreases, which 
explains their weak IRAS 60 /im emission. Calculations, for in- 
stance using the discrete dipole ap proximation (PDA) m ethod, 
of optical properties of aggregates (Bazell & Dwekl [l990) show 
that BG aggregates are more efficient emitters than a collection 
of individual grains. This increased emissivity leads to more ef- 
ficient cooling for a given ISRF intensity and therefore explains 
the lower temperature observed. The calculations indicate that 
the aggregate emissivity should progressively deviate from that 
of the individual grains in the IR (around 20 //m), and remain 



parallel to it at longer wavelengths, in particular over the entire 
submillimeter and millimeter range. 

The adjunction of VSGs in the aggregate seems to have lit- 
tle effect on FIR/mm properties and simply provides a slightly 
higher absorptivity of t he aggregate i n the near-IR dStepnikl 
1200 lab . Calculations by I Stepnik et ail d2003l) show that a total 
number of at least 20 individual BG in each aggregate is neces- 
sary to provide the required emissivity increase. However, for the 
aggregation process to be efficient, the individual grains should 
probably be covered with an ice mantle, to increase the sticking 
efficiency. 

The parameter e v in Eq.Q] is consistent with the definition 
used in iBoulanger et al.l d!996l) . The emissivity e v is related to 
the absorption efficiency Q a b s , which is directly related to the 
refractive index of the material composing the grain by 



e v (A) = na Q ahs (A, a)X d , 



(3) 



where a is the grain radius. Therefore, for a single grain size, the 
wavelength dependence of e v and Q a b s are the same. The purpose 
of this paper is to determine the spectral dependence of the dust 
absorption cross section in the atomic and molecular phases of 
the ISM and to search systematically for evidence of emissivity 
changes in the dense molecular medium, to check whether the 
results obtained for a limited number of isolated clouds can be 
generalised. In the following, we use Eq.Q] to derive the wave- 
length dependence of the dust emissivity, by simply dividing the 
observed emission l v {A) by B V {A, Td)Nn, for a single tempera- 
ture Td for each ISM phase. However, we note that, for a realistic 
grain population, Eq.Q]must be integrated over the grain size dis- 
tribution. Since the dust temperature depends on the dust grain 
size for a given radiation field, the black-body term in Eq.Q]can- 
not be factorised out from the integral. Similarly, most lines of 
sight (LOS) must exhibit variations in the radiation field, lead- 
ing to a mixture of dust temperature along the LOS. As a con- 
sequence, Eq.Q]only measures an "effective emissivity", whose 
wavelength dependence may differ from that of the true emissiv- 
ity of the material constituent of the grain. In Sect.|7] we show 
however that our approach is valid over the range of dust tem- 
perature and wavelengths considered here, and that the effec- 
tive emissivity precisely recovers the wavelength dependence of 
the true material emissivity in the presence of both a realistic 
grain size distribution and a realistic mixture of the radiation 
field along the LOS. In the following, we therefore simply refer 
to £ v as the dust emissivity, and use Eq.Q]to derive the emissivity 
directly. 

In Sect.|2] we present the data that we used for our study and 
in Sect. [3] we define the emission to column density conversion 
factors used. In Sect.|4] we describe the correlation procedure 
between dust emission and gas tracers. In Sect.|5] we explain 
how we compute the grain emissivity. The results of this study 
and the verification of our hypothesis of a single temperature 
along the LOS are presented in Sects. [6] and [7] Sections|8]and|9] 
are devoted to our discussion and conclusions. 



2. Data 

2.1. FIR data 

The Diffuse Infrared Background Experiment (DIRBE) was 
an infrared photometer onboard the COBE satellite (launched 
in 1989) to measure the diffuse infrared and microwave ra- 
diation from the early universe. It observed the entire sky at 
10 different wavelengths between 1 and 240 //m, with a 40' 
instantaneous angular resolution. However, the asymmetric 
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beam of the DIRBE instrument convolved with the spinning 
of the instrument produced an effective beam of g 1° in the 
yearly averaged products dCambresv et all |2001b). For our 
study, we only consider data in the far-infrared at 100, 140, 
and 240 yum, since emission at shorter wavelengths contains 
a large contribution from thermally fluctuating VSGs and 
polycyclic aromatic hydrocarbons (PAH). The DIRBE data 
at these wavelengths were calibrated with models for Uranus 
(100 /im) and Jupiter (140 and 240 fim) together with in-flight 
beam shape measurements. This calibration matched that of the 
FIRAS instrument, which incorporated an absolute calibrator, 
within a small correction factor, which we did not a pply in this 
study since it is smaller than our uncertainties dFixen et al.L 
11997b . 

Archeops was a balloon-borne experiment dedicated to the 
measure of the temperature fluctuations in the CMB. The focal 
plane instrument, a multi-band photometer, worked in four 
bands centered on 550 //m, 850 pun, 1.4 mm, and 2.1 mm . 
Arche ops had an angular resolution of 8' (see iBenoit et all 
2002, for a full description of the instrument). We use data 
obtained during the last flight of the instrument from Kiruna 
in February 2002, which covered about 30% of the sky. We 
note that we limit our analysis to the region surveyed by 
Archeops, which is given by the range of Galactic longitudes 
75° < 1 < 198°, in the second and third Galactic quadrants. The 
Archeops data were calibrated against the FIRAS maps, at 550 
//m and 850 fim, and with respect to the CMB dipole at 1.4 
mm and 2.1 mm. We note that the calibration of the Archeops 
data should be accurate at all wavelengths, since it relies on 
the FIRAS maps at high frequenci es, which are absolutely 
calibrated to 3% dFixsen et all ll 9941) . and on the CMB dipole 
at low frequencies, whose brightness is accurately known. In 
addition, the Archeops cosmology results derived from the 
low frequency channels have been shown to be consistent with 
those of the WMAP satellite at the 7% accuracy dTristram et all 
2005b. Details of t he A rcheops data processing are given in 



Maci as-Perez et all d2007l) . We note that the Archeops channel 



at 550 /im has residual stripes, because only one detector was 
available at that frequency. 

The Wilkinson Microwave Anisotropy Probe (WMAP) mea- 
sured the emission over the entire sky in the microwave range 
(3.2 - 13 mm), and provided accurate maps of the CMB fluctu- 
ations. The WMAP data, with a 13' beam, is currently the in- 
strument of the highest angular resolution covering the entire 
sky in this wavelength domain. Here, we use only the W band 
data (3.2 mm), because the other bands are often dominated 
by gas emission, such as free-free a nd synchrotron, or by the 
anomalous foregrou nd emission (see Draine & Lazarianl [1998; 
Benn ett et all l2003). presumably because of small spinning par- 
ticles. The calibration of the W MAP data was comp leted using 
the CMB dipole amplitude (see lHinshaw et alll2003l) . 

2.2. Gas tracers 



We us e the Leiden/Dwingeloo survey by lHartmann & Burtonl 
(1996) for the HI data, with sky observations above -30° of 



Galactic latitude, obtained with the 25 m Dwingeloo telescope, 
whose a 
details). 



whose angular resolution is 36' (see iHartman 



ingeloo t 
ni l 1994 



for more 



We use 12 CO(J = 1-0) data compiled by iDame et al] 
(12001b . These observations were obtained along the Galactic 
plane in a Galactic latitude range from 4° to 10°, with an angular 



resolution of 7.5', and some observations of large clouds at 
higher latitude, with an angular resolution of 15'. The spatial 
coverage is about 45% of the sky. 



All data were projected onto the HEALPix pixelisation 
scheme (Hierarchical Equal Area isoLatitude Pixelisation) 1 with 
a nside=128, corresponding to a pixel size of 0.45°. For the 
Archeops and WMAP first release data, we use the published 
maps, and for HI and CO, we use the maps available on the 
WMAP lambda web site 2 . For the DIRBE data, we used our 
own resampling method developed in the context of the ancil- 
lary data for the Planck mission, computing the pixel intersec- 
tion between the original DIRBE maps in the sixcube format and 
the HEALPix pixelisation. The method preserves the photome- 
try without significantly affecting the angular resolution. Since 
the Archeops data were filtered to subtract slow drifts, we sim- 
ulated fake data stream with all data described in this section, 
according to the actual Archeops scanning strategy. Those data 
time-lines were treated in a similar way to the Archeops ones, in 
particular by including the low frequency filtering applied to the 
Archeops data, and were then reprojected into sky maps. Using 
the standard HEALPix tools, all the data were corrected to the 
DIRBE angular resolution (1°), by convolution with a Gaussian 
kernel of an appropriate size in relation to the original resolution 
of each dataset. 



3. Conversion factors 

Assuming that the gas is optically thin, the hydrogen column 
density can be deduced from the integrated intensity of the HI 
emission at 21 cm (Whi) using the relation 



N m = X m W H i, 



(4) 



where Xhi is the HI integrated intensity to column den- 
sity conversion factor, assum ed to be equal to 1.82 x 

10 18 H 2 /cm 2 /(Kkm/s) (SpitziE fl978l) . The CO to molecular 

N H 

conversion factor is given by Xco = where Wco is the 
integrated CO emission, and we assume a fixed value of Xco, 
which is def i ned to be the standard Galactic value proposed by 
IStrong et al.1 d 1988b 



X co = 2.3 x 10 20 H 2 /cm 2 /(Kkm/s). 



(5) 



This value was inferred from y-ray measurements over the entire 
Galactic plane excluding the Galactic center region. 



4. Correlation between IR and gas tracers 

We considered only regions along the Galactic plane with |b| < 
30° and pixels with sufficient HI or CO emission (Whi > 
lOOOK/km/s and Wco > 0.5K/km/s). To assess the possi- 
ble contribution from free-free, synchrotron, and spinning dust 
emission on the IR emis sion in the WMAP W band, we used 
the Planck Sky Model 3 (iDelabrouille et all l2009h . The model 
indicates that the free-free emission is the main contribution in 



1 HEALPix pixelisation distributes 12 x Nside 2 points as uniformly 
as possible on a sphere surface, knowing that these points are divided 
into (4Nside-l) parallel on latitudes and are fairly spaced on longitudes 
on each one of these parallels. See http://healpix.jpl.nasa.gov/ 

2 (http://lambda.gsfc.nasa.gov/) 



see http://www.apc.univ-paris7.fr/APC_CS/Recherche/Adamis/PSM/psky.php 
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this band and generates about 30% of the IR emission. Its con- 
tribution was subtracted from the W band emission, using the 
free-free map predicted by the Planck Sky Model. Other con- 
taminations (synchrotron, spinning dust) are negligible in this 
band (< 3%), according to the same model. The contributions 
from free-free, synchrotron, and spinning dust are insignificant 
in other data used in this study. 

We performed correlations between the infrared emission 
and gas tracers, in a set of rectangular regions covering the area 
of interest. Individual regions have sizes Al x Ab= 6° x 4°, dis- 
tributed across a regular grid of Galactic coordinates every 3° in 
longitude and 2° in latitude. In each region and each map, we 
subtracted a background value computed as the median over a 
common background area, defined to be the faintest half of the 
HI data. In this way, we ensured that the correlation produces 
a null IR emission for a null column density. This step also re- 
moves any possible residual contribution from both the Zodiacal 
light and the cosmic infrared background. We then determined 
the best-fit linear correlation between the FIR emission and gas 
tracers using 



X. 



/// 



KiA) 
2X, 



N C H ° + c v (A), 



CO 



(6) 



where I v is the IR emission brightness at wavelength A, N™ and 
are the hydrogen column density in the atomic and molec- 
ular phases, respectively, and c v is a constant. The constant term 
is usually small but could account for the contribution of an ad- 
ditional gas component unaccounted for in our description, such 
as diffuse ionized gas, provided it is partially decorrelated from 
the atomic and molecular gas distribution. We performed a x 2 
minimization searching for a v and b v at each IR wavelength and 
for each region using the interactive data language (IDL) linear 
regression function "regress". 

In a second step, we rejected regions of smaller than 15 pix- 
els to maintain a sufficient number of pixels to perform the cor- 
relations. The number of pixels in each region is given in tables 
in the Appendix. With this selection, we removed 15% of the re- 
gions. In a third step, we removed regions with clearly unphys- 
ical values of the dust temperatures (> 1000 K). For tempera- 
tures inferred from the spectral shape of a v and b v between 100 
and 550 fim (see Sect.|5]for more explanation), these unphysical 
values imply that at least one of the correlation coefficients a v or 
b v in this wavelength range is still affected by uncorrected instru- 
mental effects, mostly residual stripes in the Archeops 550 fim 
channel. This selection removed 27% of the remaining regions. 

5. Grain emissivity determination 

5.1. Method 

Using Eqs.Q]and[6] an d assuming a single temperature of dust in 
the atomic (T™) and the molecular (T9°) phases, the emissivity 
of the dust associated with each phase can be written 



J'(A) 



ay(A) 1 
B v (T™)X m ' 



and 

e co ( , )= b v (A) 



B v (T c d °) 2X C0 ' 



(7) 



(8) 



respectively. The values of the dust temperatures in the two 
phases can be obtained by fitting the SED of the correlation co- 
efficients a v and b v , which is equivalent, for a given column den- 
sity, to the emission spectrum of dust associated with the atomic 



and molecular phases, respectively. We fitted the above SEDs 
by searching simultaneously for the best-fit average temperature 
and dust emissivity index for each phase, using the correlation 
coefficients at 100, 140, 240, and 550 /mi. Uncertainties in the 
emissivity values of each region are computed using 



. . AB V (T) M v 

^ = w + - 



(9) 



where the uncertainties on the intensity for the atomic and 
molecular phase are given by AI v hi = ^ and ALro = 
respectively. Errors in the a v and b v parameters correspond to 
the 1-cr standard deviation computed over all valid pixels of the 
region. Given the Planck function 



B V (J) 



2hv 3 



1 



1 



(10) 



the relative error in B V (T) (AB V (T)) is related to the error in T 
(AT) by 



AB V (T) ffes? AT 



B Y (T) 



em 



1 



(ID 



Errors in temperature used to derive the error in the Planck func- 
tion are determined using a bootstrap method, when fitting for 
the temperature and spectral index. We note that the tempera- 
tures are mostly constrained by values close to the peak of the 
emission (i.e., 100, 140, and 240 /mi). However, uncertainties in 
the derived temperature will strongly affect the emissivity values 
close the peak, because we divide by the black-body function, a 
non-linear function of the temperature near the peak. Therefore, 
the uncertainty in the temperature is the main source of uncer- 
tainty in this study. Overestimating the dust temperature by 2 K 
leads to an overestimate of B v and therefore an underestimate of 
the dust emissivity by as much as 150% and 70% at 100 fim and 
240 /an, respectively, for an average dust temperature of 16 K. 
For AT = 0.5 K, these figures become 65% and 35%, respec- 
tively. 

All data used in this study follow a given flux convention 
(e.g., vL = constant for DIRBE data), which allows us to com- 
pute unambiguously the total power received in the instrument 
photometric band, regardless of the energy distribution in the 
band. However, the true brightness at the reference frequency 
depends on the emission spectral shape, which is unknown un- 
til a fit using a model has converged. We therefore corrected the 
emissivity derived above by dividing by the appropriate color 
correction factor, computed using the filter transmission and flux 
convention of each instrument used. This was done in an itera- 
tive way, starting from color correction factors derived from a 
dust model, and iterating for the true shape of each SED, until 
convergence was reached. 

5.2. Classification according to dust temperature 

Dust in dense molecular clouds is expected to be colder than 
in the surrounding atomic material, as long as no significant 
star formation is occurring in the cloud. However, massive star 
formation can significantly heat the dust in molecular clouds. 
In addition, star formation may increase significantly the tur- 
bulence in the star-forming molecular medium, which in turn, 
could modify the dust emission properties. To search for these 
potential variations, we classify the regions according to 3 cat- 
egories based upon the dust temperature derived for dust in the 
atomic phase ( T?) and in the molecular phase (T9°): 
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- regular: atomic and molecular medium with similar dust 
temperatures defined as < 1 < + 0.6 K. This case 
allows us to compare the dust emissivity in the two environ- 
ments in regions where the ISRF intensity is similar in both 
phases. This category accounts for 18 regions in our sample. 

- colder CO: dust in the molecular phase significantly colder 
than in the atomic phase: T™ + 2.5K< T™ < T^° + 4K. 
In this case, the colder dust temperature in the molecular 
phase is indicative of a lower ISRF intensity in the cloud, 
most likely indicating dense molecular environment with lit- 
tle star formation. This case is the most relevant to the search 
for emissivity variations linked to the formation of dust ag- 
gregates, sinc^_sofjrr^hesehave_b^ regions 
(see e.g jBernard et al.Lll999trStepnik et alll2003l) . This cat- 
egory accounts for 45 regions in our sample. 

- warmer CO: dust in the atomic phase significantly colder 
than in the molecular phase: T™ + 2 K < T£° < T™ + 4 K. 
This category should incorporate most of the massive star- 
forming regions. The lower temperature (here T™) is lower 
than in the previous case 2 in order to include enough re- 
gions. This category accounts for 25 regions in our sample. 
FigureQ] shows the location of all studied regions for the 
three temperature cases. Temperature histograms are shown 
in Fig. |2j for each case and phase. Although the cases are 
selected based on the relative temperature between the two 
phases, it can be seen that they also select absolute tempera- 
ture behaviors, since the absolute values of the dust tempera- 
tures are roughly reversed in case 2 and 3 (T^j Se2 T^Q Se3 =2 
17 K and Tgf 2 =* T^ se3 =* 14 K). Appendix (Table 7 and 
8) indicates for each case and studied region, the number of 
pixels considered, the correlation factors obtained as well as 
the dust temperature values, and the emissivity spectral index 
in both the FIR and the submm range. 

6. Results 

Figure|3] shows the SEDs of the median emissivity derived 
from the correlation values, for the atomic and molecular 
phases, for each of the cases described in Sect. |5.2| The cor- 
responding values are also given in Table Q] 
The emissivity values shown in Fig. [3] were not normalised 
in any other way than assuming the fixed conversion factors 
described in Sect. [3] In Figure|3] we also sh ow the emissiv- 
ity val ue derived from the FIRAS data by Boulanger et al. 
(1996) for the atomic diffuse ISM in the solar neighbour- 
hood at A = 250 yum, assuming ft = 2. It can be seen that the 
emissivity values derived here are in rough agreement with 
that value. We note however that the atomic emissivities that 
we derived are somewhat higher than in the solar neighbour- 
hood. We also note that there is a continuity in the emissivity 
SED between the WMAP W band and the Archeops range 
indicating that there is no obvious calibration difference be- 
tween the two experiments, which are both calibrated to the 
CMB dipole in this range. 

The scaling of emissivity values for the molecular phase 
can be modified by assuming a different Xco factor. The Xco 
factor is difficult to determine. It is expected to vary with var- 
ious factors such as the metallicity and the chemistry of the 
ISM, and could therefor e change from cloud to cloud (e.g., 
Mag nani & Onellolll995h . In Figure[4] we show the same emis- 
sivity SEDs, where we have arbitrarily chosen the Xco value so 
that the atomic and molecular emissivities match at the longest 
wavelengths. For each case, the derived Xco values are given in 




Fig. 1. Studied regions for each temperature case: case 1 (regu- 
lar) in white, case 2 (colder CO) in grey and case 3 (warmer CO) 
in black. The sky region covered by the Archeops data is also de- 
lineated. The coordinate grid shown is in Galactic coordinates. 
The map shows half the sky and is centered roughly towards the 
Galactic anticenter. 



the plot. Figure|4]more clearly shows the variations in the spec- 
tral shape. From that figure it is apparent that in case 1 and 3, the 
atomic and molecular medium emissivity SEDs are almost par- 
allel, while in case 2 (colder CO), the molecular emissivity SED 
is significantly steeper than the atomic one in the FIR range. We 
note however, that if the absolute level of the emissivities in the 
molecular phase can vary according to the assumed Xco value, 
the spectral shape of the emissivity SED and therefore the emis- 
sivity slope (J3) will remain unchanged. 

In Figure|U we also show power laws with /3 - 2.5, (3 — 2 
and p = 1.5 normalised at 100 /.im for the atomic phase. 
Although all emissivity SEDs globally show the expected de- 
crease with wavelength, it can be seen that, in all 3 cases, 
the emissivity SED is steeper in the FIR (DIRBE data) and 
flattens significantly in the submillimeter above A > 500 jum 
(Archeops and WMAP data). Tables [2] and [3] summarise the me- 
dian values of the emissivity spectral index at all wavelengths 
in the DIRBE and Archeops wavelength ranges, respectively. 
These values were inferred from a linear adjustment of the emis- 
sivity SEDs. Uncertainties correspond to 1-cr. It can be seen that 
the emissivity index is generally higher than ff 11 - 2 in the FIR 
corresponding to a steep emissivity decrease with wavelength, 
and is about yg submm ~ 1 .5 in the submm and mm domain, corre- 
sponding to a flatter spectrum. Apart from this common behav- 
ior, there are also significant differences between the 3 tempera- 
ture cases, which are discussed below. 

6.1. Case 1: regular 

In this case, the emissivity SED of the dust associated with 
the atomic and molecular phases are almost parallel, and both 
clearly experience the same slope change between the DIRBE 
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Table 1. Median values of the dust emissivity in the molecular and the atomic phase for each studied case: case 1 (regular), case 2 
(colder CO), case 3 (warmer CO). Error bars represent the + 1-cr uncertainty estimates. 
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Table 2. Median values of the emissivity spectral index (J3) in the molecular and the atomic phase for each studied case and each 
wavelength. Error bars represent the 1-cr uncertainty estimates. 
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Table 3. Median values of the emissivity spectral index (J3) for the FIR and submm domain, for each studied cases. Errors bars 
represent the 1-cr uncertainty estimates. 



and the Archeops wavelengths. This similarity probably indi- 
cates that, in this case, the dust emission properties associated 
with the molecular and atomic phases are the same. 

In addition, the absolute emissivity values in the two phases 
are relatively close. This indicates that our adopted Xco value 
is reasonable in that case. There is no reason in principle that 
the dust abundance Xd could differ between the two phases. 
Therefore, an incorrect estimate of the Xco parameter would 
produce a wavelength-independent difference between the two 
emissivity SEDs. Looking at the median values of the emissiv- 
ity ratios in the two phases presented in Table|4] we can see that 
the emissivity in the atomic phase is higher that in the molecular 
phases at all wavelengths, by about 30% for the assumed Xco- 
This difference could be accounted for by a change in the Xco 
value, which would then have to be slightly lower than the one 
assumed, X co = 1.6 X 10 20 H 2 /cm 2 /(Kkm/s). 

6.2. Case 2: colder CO 

This case includes regions where dust in the molecular phase 
is significantly colder than in the surrounding atomic medium. 
It is of the greatest interest when comparing our results with 
those obtained using the PRONAOS data in the Taurus and the 



Polaris regions, which lead previous authors to propose the pres- 
ence of fractal dust aggregates in the molecular phase. Our re- 
sults (see Fig. [3) indicate that the emissivity SED of dust asso- 
ciated with the atomic and molecular medium differ. In the FIR 
(DIRBE wavelength range), the emissivity spectral index in the 
molecular phase is significantly higher than in the atomic phase. 
According to Table[3]the molecular to atomic emissivity ratio 
(tt 2 ) is equal to 1.3 in the FIR domain. However, in the submil- 

Phi 1 

limeter (Archeops wavelength range), the emissivity spectra in 
the two phases become parallel again, and have quite compara- 
ble emissivity values for the assumed value of Xco- Matching 
exactly the submm emissivities in the two phases in the submil- 
limeter range would lead to X C0 = 1.5 X 10 20 H 2 /cm 2 /(Kkm/s). 
According to Table|3] ^=1.0, indicating that the slope of the 
submm emissivity is not very different from case 1, whose ratio 
is 0.9. For the adopted Xco, we observe that the emissivity in 
the molecular phase is higher than that in the atomic phase in 
the FIR range, by a factor ranging from 3.1 at 100 //m, to about 
1.5 at 240 fim (see Table |4j. The absolute emissivity values 
for the two phases tend to become progressively similar around 
550 yum. This behavior cannot be attributed to the assumed Xco 
value, which would affect emissivity values at all wavelengths. 
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Fig. 2. Dust temperature histograms for each selected case: the 
atomic and molecular phases are shown in dark grey with black 
contours and light grey, respectively. 



mo - , .lore 

Fig. 3. Median dust emissivity SEDs for case 1 (regular: upper 
panel), case 2 (colder CO: middle panel) and case 3 (warmer 
CO: lower panel). The DIRBE, Archeops and WMAP data cor- 
respond to diamond, star, and square symbols, respectively. The 
shaded areas show the ± 1-cr dispersion around each emissiv- 
ity SED, in dark grey for the atomic phase and in light grey for 
the molecular phase. The triangle symb ol at A = 250 /jm corre- 
sponds to the emissivity value derived by Boulan ger et al.l (fl996) 
for the diffuse atomic medium at 250 fim. 



We also considered the possibility that this could be due to an 
error on the determination of the dust temperature in the molec- 
ular phase. Obtaining the same emissivity slope in the FIR in 
the two phases would require an increase of the temperature in 
the molecular phase of 2-3 K, which is larger than the temper- 
ature uncertainty, which are between 0.25 K and 0.72 K. In our 
analysis, we have so far used the major simplification of a single 
dust temperature along the line of sight. In Sect. [7] we investi- 
gate the effect of this hypothesis on the results. In particular, we 
consider the effect of the temperature mixing due to the grain 
size distribution, the LOS mixing of the radiation field and the 
grain composition. 



In this case, the dust temperature in the molecular phase is 
higher than the dust temperature in the atomic phase. Taking 
into account the error bars, the behavior of the emissivity spec- 
trum is qualitatively the same as in case 1. For the assumed 
Xco, we note that the atomic emissivities are significantly 
higher than the molecular ones at all wavelengths. To repro- 
duce the emissivity in the submm, we would require Xco=l.l x 
10 20 H2/cm 2 /(Kkm/s). The results of this case show that when 
the dust temperature in the atomic phase is close to or colder 
than the molecular phase, grains emission properties seem to be 
the same. 
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0.7 ± 0.2 


0.6 + 0.2 


850 


0.6 ± 0.2 


0.5 ± 0.2 


0.5 ± 0.2 


1400 


0.8 ± 0.2 


0.8 ± 0.2 


0.5 ±0.1 


2100 


0.7 + 0.2 


0.6 ± 0.3 


0.3 ± 0.2 


3200 


0.7 ± 0.2 


0.7 ± 0.4 


0.3 ± 0.3 



Table 4. Median values of the ratio between the dust emissivity 
in the molecular phase and in the atomic phase for each stud- 
ied case. Error bars represent the + l-<x dispersion around each 
value. 



7. Possible biases caused by assumptions 

In the study presented above, we assumed a single dust tempera- 
ture along each LOS. This is probably overly simplistic, since we 
expect a mixture of temperatures for various reasons. First, dust 
exhibits a size distribution, and dust temperatures depend upon 
grain size for a given radiation field. Second, the ISRF strength 
probably varies along any given LOS within the Galactic plane, 
even towards the external Galactic regions studied here. Third, 
composition variations in the dust could produce temperature 
variations, which could also affect our results. Since we divide 
the observed sky brightness by the black-body function at a sin- 
gle temperature, slight systematic biases in the dust temperature 
caused by these effects could in principle explain the inferred 
variations of the dust emissivity, particularly in the FIR region, 
where the black-body function is non-linear with respect to the 
dust temperature. 

In this section, we explore those 3 possibilities and test the 
robustness of our results to temperature variations along the LOS 
induced by either grain size distribution, ISRF strength mixing 
or grain composition variations. In all cases, we produce pre- 
dicted emission SEDs using a pertinent dust emission model, 
which includes the additional complexity, apply the same treat- 
ment applied to the true data to the modelled SEDs, and compare 
the derived emissivity curves with that of the model. 



7. 1 . Case 3: warmer CO 

7.2. Effect of the grain size distribution 

To test the influence of temperature caused by the BG size dis- 
tribution, we generated model SE Ds at various ISRF intensities 
(Xtsrf), using the dust model of Dese rt et"aT] 0990), which in- 
cludes a realistic size distribution for large grains and computes 
the induced temperature mixing, by taking into account the ef- 
fects of grain temperature fluctuations. Model SEDs were com- 
puted in all instrument filters used in our study, using the proper 
color correction to match the flux convention of each instrument. 
We then applied the same processing to the model SEDs as ap- 
plied to the actual data. As the model assumes a dust emissivity 
with /3 — 2 over the entire wavelength range considered here, 
we expect in that case to recover a simple power law with slope 
P = 2. We performed the test over Xisrf values ranging from 
0.05 to 10 times the solar neighbourhood ISRF strength. Ratios 
of the predicted emissivity values inferred from the model and 
the apparent emissivity obtained with our procedure are given in 
Table[5] normalised to the emissivity value in the WMAP band. 




mo - , .lore 

Fig. 4. Same as Fig.[3]but with the molecular emissivity scaled 
to match that of the atomic phase in the range 550 fim < A < 
3 mm. The corresponding Xco values are given in the figure. For 
comparison, power law spectra with /J = 1.5 (dashed), = 2.0 
(solid) and (3 — 2.5 (dot) are shown, normalised to the atomic 
emissivity at 100 fim. 



Figure [5] is a graphical representation of Table [5] The left panel 
shows that the emissivity ratio is close to 1 . We therefore recover 
emissivity spectra very similar to the input power law, indicating 
that no significant bias is introduced by not considering the LOS 
temperature mixing induced by a realistic dust size distribution. 
We also checked the influence of changing the spectral shape of 
the ISRF, in particular increasing the UV content as expected 
close to young stars. This test showed that the shape of the ISRF 
does not affect our results, because BGs absorb uniformly across 
the entire visible to ultra-violet range of the model used. We con- 
clude that the method described above introduces no significant 
bias in that case. 



D. Paradis et al.: Dust emissivity variations in the Milky Way 



9 



A (fJ.m) 


Grain size distribution 




ISRF strength mixture 


Grain composition 






X IS rf=0.05 


Xisrf- 0.5 X 1SRF -4 


a=2.5 


a= 2 


a= 1.25 


12 = 1UK 12= 15 K 1 


_=Ll IS. 


100 


1.13 


0.89 


0.93 


1.20 


1.46 


3.06 


1.06 1.29 


1.16 


1 Af) 


1.13 


0.96 


0.98 


1.23 


1.45 


2.66 


11? 1 77 


1 14 


240 


1.11 


1.00 


1.02 


1.17 


1.29 


1.84 


1.03 1.19 


1.13 


550 


1.07 


1.03 


1.03 


1.09 


1.12 


1.27 


0.96 1.05 


1.04 


850 


1.04 


1.02 


1.02 


1.05 


1.07 


1.14 


1.07 1.10 


1.08 


1400 


1.00 


0.99 


0.99 


1.01 


1.01 


1.05 


1.00 1.04 


1.05 


2100 


0.99 


0.99 


0.99 


0.99 


1.00 


1.01 


0.98 1.01 


1.03 



Table 5. Ratios of the predic ted to apparent emissivities, using the same processing as for the data with different models. Cols 2-4: 
using the Dese rt et ail d 19901) model at various IS RF strength (Xisrf). Cols 5-7: assuming a mixture of the radiation field intensity on 



each LO S, using the model bvlDale et a for various values of the mixing parameter a. Cols 8-10: using the two component 



model of Finkbeiner et alj d 19991) for different temperatures of the warm component (T2). 
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Fig. 5. Ratios of the predict ed (e p ) to apparent emissivities (e), using the same processing as for the data with different modelled 
SEDs. Left panel: using the lDesert et al.l (1 19901) model at various ISRF strengths (Xisr f)- Middle panel: a ssuming a mixture of the 
radiation field intensity on each line of sight, using the model descr iption proposed bvlDale et al.l d200ll) for various values of the 
mixing parameter a. Right panel: using the two component model of Finkbei ner et al.l d 19991) for different temperatures of the warm 
component (T2). 



7.3. Effect of the ISRF strength mixture 

To test the influence of the variations in the ISRF str ength mix- 
ture al ong the LOS, we use the model proposed by iDale et alj 
(2001), who introduced the concept of local SED combination, 
assuming a power-law distribution of dust mass subjected to a 
given heating intensity dMd(XisRF), to interpret the emission by 
external galaxies: 



dM d (X 1SRF ) °c X"dX 1SRF , 



(12) 



where a controls the relative co ntribution o f the v arious ISRF 
intensities (Xisrf) to the SEDs. IDale et al.l d200ll) used ISRF 
strengths values in the range 0.3 < Xisrf < 10 5 and derived 
a values in the range 1 < a < 2.5 for external galaxies. W e first 
computed emission spectra using the lDesert et al.l dl990h model 
at various Xisrf values (noted I™° d (XisRF )), then summed them 
over the same Xisrf range as proposed bv lDale et al.l d200ll) . ac- 
cording to 



jtot _ 



ISRF 



1 Y~ a 
■iXjsgF ISRF 



(13) 



In Figure|5] (middle plot), it can be seen that significant depar- 
tures from model predictions can occur in the FIR, close to the 
peak of the emission, for flat Xisrf distributions (i.e., low a 
values). These deviations become progressively smaller towards 
steeper distribution, as expected. We note that the apparent dust 



temperature derived from the fit decreases as a increases, which 
can be used to constrain a values. From the model predictions, 
we corrected our emissivity estimates, using the apparent dust 
temperature for each measured SED to infer the true a value 
and then deduce the factor we have to apply to correct our val- 
ues, knowing the ratio of the predicted to apparent emissivity. 
In general, assuming a single dust temperature along the LOS 
causes the emissivity values between 100 and 550 fim to be un- 
derestimated. For most of the regions considered here, the appar- 
ent temperature is moderate, which is indicative of steep Xisrf 
mixing with a close to 2.5. The resulting corrected median emis- 
sivity values are shown in Fig. [6] after renormalisation by Xco- 
Comparison with Fig. [4] shows that Xisrf mixing has little im- 
pact on our results. In particular, both the break in the emissivity 
slope near 500 /vm and the dust emissivity increase in the molec- 
ular phase in case 2 are still present. We note that the amplitude 
of the emissivity increase is slightly lower after the correction, as 
shown in Table|6] with a median ratio of 2.0 ± 1.6 between the 2 
spectra at 100 //m after correction, instead of 3.1, before renor- 
malisation by Xco- The emissivity slopes in the FIR are equal to 
2.7 + 0.6 and 2.3 + 0.6, in the FIR for the molecular and atomic 
phase, respectively, and 1.5+ 0.2 and 1.6 + 0.2 in the submm. 



7.4. Effect of the grain composition 



We use the Finkbei ner et al.l (|T999) model to test the influence of 
dust composition on our results. In the framework of this model, 
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Table 6. Median values of the dust emissivity in the molecular and the atomic phase for each studied case, after correction of a 
mixture of the radiation field intensity along the line of sight. Error bars represent the + 1-cr uncertainty estimates. 



the shape of the submm emission is assumed to be caused by a 
mixture of silicate and graphite grains that reach different tem- 
peratures (Ti and T2 for the silicate and graphite components, re- 
spectively), the temperatures being linked by the radiation field. 
In principle, this could also affect our interpretation, since we 
assumed a single du st temperature to deriv e the dust emissiv- 
ity spectrum. Unlike Finkbeineret al. (1999) however, we used 
the same A~ 2 power law emissivity for each grain component, to 
simplify the comparison between the recovered emissivity and 
the model. We computed the emissivity ratio at various temper- 
atures of the graphitic component ranging from 8 K to 22 K (see 
the right plot of Fig.[5j- Table[5]shows that in this case the bias in- 
troduced by our simplification hypothesis does not significantly 
impact our results. Similar tests using various emissivity power 
laws for the two components gave equally good results. 

8. Discussion 

Regarding the absolute values of the derived emissivities, it can 
be seen in Fig.[3]that the values we derived for t he atomic phase 
are in rough agreement with that derived by Boulanger et al.l 
dl996l) (1(T 25 cm 2 /H at 250 /urn). This agreement indicates that 
the regions we have selected in this study, which are mainly lo- 
cated towards the Galactic 2 nd quadrant, have similar dust con- 
tent and properties to the solar neighbourhood. A closer inspec- 
tion of Fig. [3] shows that our emissivity values in the atomic 
phase are systematically higher than the solar neighbourhood 
values by a factor of 1.9, 1.3, and 2.6 for cases 1, 2, and 3, re- 
spectively. The corresponding median dust temperatures derived 
are 15.0 K, 16.9 K, and 14.0 K respectively (se e Table 7), which 
are als o lower than the temperature inferred by Boulanger et alj 
(1996) for the diffuse ISM in the solar neighbourhood (17.5 K). 

We also emphasize that the atomic and molecular emissiv- 
ity spectra for cases 1 and 3 could be reproduced over the full 
wavelength range studied by slightly decreasing the Xco value 
used for the column density conversion. We therefore consider 
it likely that the actual Xco value is slightly lower than the as- 
sumed value of 2.3 X 10 20 H2/cm 2 /(Kkm/s) in those cases. In 
case 2, however, the emissivity spectra have a significantly dif- 
ferent shape and cannot be reproduced over the full wavelength 
range by simply changing the Xco value. We note however that 
a similar decrease to that in case 1 would improve the match in 
the submm and millimeter range. 

Regarding the spectral shape of the emissivity, the results 
shown in Sect. |6. 21 indicate that the FIR dust emissivity of dust 
associated with molecular regions is significantly higher than 
that of the surrounding atomic medium, in cases where the dust 
temperature is colder in the molecular phase than in the atomic 



phase (Case 2). Taking into account the error bars derived in 
Sect. l5.1l and shown in FigHl this result is significant at the 5 - <x 
level. Thi s result is in qualitat ive agreement with the result ob- 
tained by Step nik et al.l d2003l) . who showed that an increase in 
the dust emissivity by a factor of 3-4 in the molecular phase is re- 
quired to explain the low dust temperature in that phase. We note 
that their study was limited to the wavelength range from 200 to 
600 fim. The increase that we found for the DIRBE wavelength 
range is therefore of the same ord e r, alth ough slightly smaller, 
than that derived by IStepnik et alj d2003l) for a particular cold 
molecular filament in Taurus. In our case 2 regions, we note that 
on average the dust is 3.2 K colder in the molecular phase than 
in the atomic phase. This is also simi lar to, although sligh tly 
lower than, the difference observed by IStepnik et alj d2003l) in 
the Taurus filament (4.8 K). Our result ther efore confirms the 
hypot hesis proposed by Ber nard et alj d 19991) and Stepni k~et al.l 
d2003l) that dust aggregation may be responsible for the de- 
crease in the dust equilibrium temperature. In quiescent molec- 
ular clouds of case 2, aggregation of BGs in fractal clusters is 
a likely cause of the FIR emissivity increase, the correspond- 
ing decrease in the dust equilibrium temperature being a conse- 
quence of the higher emissivity and more efficient cooling. 

Stepn ik et alj |2003) inferred an increase in emissivity by a 
constant factor with wavelength. Their data did not allow them 
to measure directly the spectral shape of the emissivity, and they 
used an average /3 index over the FIR wavelength range of the 
PRONAOS experiment. Their approximation of a wavelength in- 
dependent increase in the emissivity was supported by the results 
of dis crete dipole appr oximation calculations for dust aggregates 
(e.g., iBazell & Dwekl[l990l) that predict a constant absorptivity 
increase for fractal aggregates all the way into the millimeter. 
Our results however show that the emissivity spectrum in the 
molecular phase follows the same power law as in the atomic 
phase. Although the dust emissivity in the cold molecular phase 
is significantly higher in the FIR, it recovers the atomic phase 
values in the submm and mm ranges. 

We note however that reproducing the temperature difference 
by increasing the dust emissivity only requires the emissivity to 
be higher in the FIR, where large grains emit most of their en- 
ergy. Therefore, the close agreement between the emissivities 
derived for the submm and mm does not therefore the thermal 
argument upon which the conclusion of dust aggregates is based. 
A simple calculation of the thermal equilibrium shows that the 
measured emissivity increase in case 2 would produce a temper- 
ature decrease of about 2 K, assuming that the increase below 
A — 100 fim equals that at 100 fim. This is slightly lower than 
the observed temperature decrease, which could indicate that the 
emissivity continues increasing below 100 /mi. 
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Fig. 6. Median dust emissivity SEDs for the 3 cases (see Fig. [3] 
for the description of the curves) corrected for the hypothesis of 
a mixture of different interstellar radiation field along the line of 
sight. The molecular emissivity has been scaled to match that of 
the atomic phase in the range 550 /mi < A < 3 mm. The triangle 
corres ponds to the emissivity value derived by iBoulanger et al.l 
(1996) in the diffuse medium. 



In all cases studied, the emissivity spectra derived are steeper 
in the FIR, with spectral index of /3 2.4 and flatten in the 
submm and mm regime, where the spectral index reaches f3 
1.5. Taking into account, the error bars derived in Sect. 15. H and 
shown in Fig.|4l this result is significant at the 10 - <x and 15 - cr 
level for cases 1 and 2, respectively. We note that attributing the 
flattening to calibration errors would require the calibration of 
both the Archeops and WMAP experiments to be wrong by more 
than a factor of 2 in the same direction, which is very unlikely. 
This behavior could in principle be due to the presence of two 
dust components at different equilibrium temperatures, as pro- 
posed in phen omenological models w ith two temperature com- 
ponents (e.g., Finkbeine r et all [1999). Such models have been 
invoked to explain the flatness of submm emission spectrum, but 



require the existence of a cold component with Tj - 9 K, whose 
origin in the diffuse ISM remains to be understood. However, we 
note that the spectral indices that we obtained in this work are 
in goo d agreement with those proposed in the Finkb einer et al.l 
(119991) model. 

A more physical approach consists of interpreting the 
submm flattening by in trinsic dust proper ties. To explain the 
FIR/mm dust emission, Me nv et ail (|2007) proposed a model 
based on some specific properties of the amorphous state. First, 
they considered a temperature-independent emission caused by 
excitation of acoustic lattice vibrations, coming from the cou- 
pling between the electromagnetic fields and a disordered charge 
distribution (DCD) that characterizes the amorphous nature of 
dust. Their description also takes into account resonant absorp- 
tion and relaxation processes associated with localized asym- 
metric two-level systems (TLS) in the grains, which produce ad- 
ditional emission with emissivity that depends on both tempera- 
ture and wavelength. Within this model, the emissivity spectral 
index is therefore predicted to change with both wavelength and 
temperature. The behaviour of a flattening of the spectrum in the 
millimeter range is in qualitative agreement with their model, 
although the data do not allow us to test further a possible de- 
pendence with dust temperature, since the temperature range 
that we can sample here is limited. Finally, we note that in this 
model, the TLS phenomenon responsible for the emission above 
A =a 600 /im is of a different nature than the vibration of DCD in 
the amorphous lattice producing the FIR emission. It is therefore 
possible that fractal aggregates consisting of amorphous individ- 
ual grains also exhibit a change in properties between the FIR 
and the submm. In particular, the TLS phenomenon operating at 
atomic scale may be less sensitive to the dust grains being gath- 
ered into aggregates, while the DCD vibration, which is a global 
phenomenon, is expected to exhibit an excess emission, as pre- 
dicted by classical DDA calculations. This may offer a possible 
physical interpretation to the emissivity increase of dust in the 
molecular phase as being limited to the FIR region. Obviously, 
detailed modelling of the interaction between the electromag- 
netic wave and a fractal aggregate consisting of amorphous ma- 
terial is needed to further investigate this issue. 

9. Conclusions 

We have analysed the dust emission from the outer Galactic 
plane region using DIRBE, Archeops, and WMAP data from 
100 //m to 3.2 mm. We have performed a correlation study of 
the FIR-mm emission with gas tracers in individual regions, and 
derived the average equilibrium temperature of large dust grains 
in both molecular and atomic phases in a set of regions along the 
Galactic plane. We used this temperature to derive the emissivity 
spectra for each phase and region. We classified regions into 3 
classes, according to the relative temperature of the dust associ- 
ated with the molecular and atomic phases. For each class, we 
derived the median emissivity profile. 

We showed that the emissivity spectra are always steeper 
in the FIR (A < 600 /mi) and flatten in the submm and mm. 
In regions where dust is significantly colder in the molecular 
phase than in the surrounding atomic medium, we produced 
an increase in the emissivity by a factor of 3 in the FIR. 
However, we showed that the emissivity increase is restricted to 
the FIR range; the emissivity spectra for the dust in the atomic 
and molecular phases become comparable again in the submm 
and mm wavelength range. 

The observed break in the emissivity spectrum, which ap- 
pears to be a general phenomenon, is consistent with the phe- 
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nomenological all-sky fit of the FIRAS data of Finkb elner et al.l 
( 1999(). It is also in qualitativ e agreement with the dust emission 
model of iMenv et alj d2007l) . which invokes quantum effects in 
amorphous solids to explain the flatness of the observed submm 
emission spectrum and also produces a break in the emissivity 
slope around 600 jjm. 

We interpret the FIR emissivity increase in molecular clouds 
containing cold dust as being caused by the coagulation of large 
grains into fractal aggregates. Previous ob servations, obtaine d 
over much smaller portions of the sky bv lStepnik et alj (|2003), 
showed that dust aggregation could lead to an increase in the 
dust emissivity by a factor of 3-4, to explain the unusually low 
dust temperatures observed. The fact that emissivities do rec- 
oncile in the submillimeter is not in agreement with DDA cal- 
culations of simple aggregates and the physical reason for this 
remains unclear. We however propose that this may be related to 
the amorphous nature of the grains contained in the aggregates. 

Finally, the absence of a detectable FIR emissivity increase 
in regions where the dust temperature in the molecular phase is 
similar or higher than that in the atomic phase is interpreted as 
the absence of fractal grains in those environments. This is prob- 
ably because the formation of such aggregates is prevented by 
the star formation activity, while it is favoured in more quiescent 
regions where turbulence is milder. 
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Table 7. Results of the IR/HI correlations in the atomic phase for all sky regions. Cols 1-2: central position of each studied regions in degree. Col 3: number of pixels used for 
the correlations in each studied regions. Cols 4-11: results of the correlations in MJy/sr for 10 21 H/cm 2 . Cols 12-13: emissivity spectral indices in the range 100 - 240 //m (f3 F1R ) 
and in the range 550 - 2100 jjm (yj submm ). Col 14: derived dust temperature in Kelvin. Errors bars represent the l-cr uncertainty estimates. 
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15.3 ± 0.5 


148.23 


-6.34 


38 


2.6 ± 0.5 


10.3 ± 1.6 


1 1.8 ± 1.2 


0.5 ± 0.2 


6.4 ±1.0 


20.0 ±1.7 


6.6 ± 


0.8 


8.0 ± 2.9 


4.8 


± 


0.6 


0.9 ± 0.2 


10.3 ± 0.2 


148.23 


2.01 


111 


3.0 ± 0.3 


9.6 ± 0.7 


9.8 ± 0.5 


1.6 ± 0.1 


5.4 ± 0.4 


11.7 ± 0.6 


3.1 ± 


0.3 


6.7 ± 0.8 


2.4 


± 


0.5 


1.6 ±0.1 


14.1 ± 0.4 


151.29 


-0.07 


110 


5.3 ± 0.8 


10.9 ± 1.3 


9.6 ± 0.8 


1.4 ± 0.1 


3.0 ± 0.4 


9.1 ± 0.7 


2.4 ± 


0.4 


6.5 ±1.2 


2.1 


± 


0.5 


1.6 ± 0.2 


17.3 ± 0.6 


i j i .zy 




1 MM 
1UU 


O 1 -i- M 1 
Z. 1 ± U. 1 


J.y ± U.j 


J. 1 ± u.z 


1 -i- M 1 

1 .z ± U. 1 


O A _i_ M O 
Z.4 ± U.Z 


/.J ± U.Z 


1 . 1 ± 


M 1 


/.4 ± U.O 


O M 

z.u 


+ 


M A 
U.4 


1 Q _i_ A 1 

1.0 ± U. 1 


1 A M -i- M A 
10. U ± U.O 


1 <A 16. 

154.36 


1 A <0 

-10.52 


51 


i o i a 1 
1.2 ± 0.1 


3.9 ± 0.2 


A A i A O 

4.4 ± 0.2 


A A i A A 

0.6 ± 0.0 


O 1 i A 1 

Z.3 ±0.1 


O 1 i A O 

o.l ± 0.3 


2.2 ± 


A 1 
0.1 


H 1 i A C 

1 .5 ± 0.6 


O A 

2.4 


± 


A A 
0.4 


11 i A 1 
1.1 ±0.1 


1 1 O i A A 

13.0 ± 0.4 


157.43 


12.46 


33 


2.4 ± 0.4 


7.7 ± 1.5 


8.4 ±1.2 


1.4 ± 0.4 


3.2 ± 0.8 


6.5 ±1.4 


0.4 ± 


0.9 


<13.3 


2.3 


± 


0.6 


2.1 ± 0.5 


13.9 ± 0.4 


160.50 


-8.43 


112 


1.4 + 0.2 


4.5 ± 0.5 


5.8 ±0.5 


0.6 ±0.1 


2.8 ± 0.3 


7.3 ± 0.6 


1.7 ± 


0.3 


6.9 ± 1.2 


2.9 


± 


0.5 


1.3 ±0.2 


12.5 ±0.3 


160.50 


2.01 


117 


2.6 ±0.2 


7.0 ±0.5 


7.8 ±0.5 


1.1 ±0.1 


3.1 ±0.5 


7.7 ± 0.8 


2.1 ± 


0.3 


5.1 ± 1.6 


2.1 


± 


0.4 


1.6 ±0.2 


14.9 ± 0.4 


163.57 


8.28 


42 


1.1 ±0.3 


2.9 ±0.9 


4.2 ± 0.8 


0.3 ±0.1 


< 1.2 


<3.1 


<0.2 






2.8 


± 


0.7 




12.7 ±0.3 


163.57 


10.37 


59 


0.9 ±0.3 


2.9 ± 0.9 


3.5 ±0.7 


0.2 ±0.1 


<1.4 


< 3.3 


< 0.4 






3.8 


± 


0.7 




11.4 ±0.3 


166.63 


8.28 


34 


1.4 ±0.3 


3.8 ± 1.0 


4.6 ± 0.9 


0.4 ±0.1 


0.3 ± 0.5 


1.2 ± 1.2 


<1.2 






2.6 


± 


0.7 




13.3 ±0.3 


172.77 


-10.52 


110 


1.3 ±0.1 


4.5 ± 0.4 


5.6 ± 0.4 


0.7 ±0.1 


2.6±0.3 


7.8 ± 0.4 


1.6 ± 


0.2 


5.9 ± 0.9 


2.7 


± 


0.5 


1.4 ±0.2 


12.7 ±0.3 


Median 


values 






































2.6 ±0.3 


7.7 ± 1.0 


7.8 ±0.8 


1.2 ±0.1 


3.1 ±0.5 


8.0 ±0.8 


2.3 ± 0.4 


6.9 ± 1.6 


2.4 ± 0.6 


1.6 ±0.2 


14.9 ± 0.4 


Case 2 


71.53 


4.10 


64 


29.4 ± 2.3 


57.5 ±4.7 


34.9 ± 3.2 


4.8 ± 0.3 


< 1.6 


24.7 ± 1.5 


7.1 ± 


0.8 


0.3 ± 3.2 


2.6 


± 


0.5 


1.8 ±0.2 


17.7 ±0.7 


74.59 


2.01 


86 


12.4 ± 3.2 


29.1 ±4.0 


22.0 ± 1.7 


2.6 ± 0.2 


4.8 ± 0.8 


18.2 ±0.9 


4.8 ± 


0.4 


9.4 ± 1.1 


2.6 


± 


0.6 


1.6 ±0.1 


15.4 ±0.5 


77.66 


4.10 


104 


41.6 ±3.3 


63.4 ±4.0 


34.8 ± 1.7 


2.6 ± 0.2 


11.9 ±0.4 


24.8 ± 0.8 


8.6 ± 


0.4 


7.3 ± 1.4 


2.9 


± 


0.4 


1.2 ±0.1 


17.5 ±0.6 


80.73 


8.28 


56 


28.2 ±5.0 


56.6 ± 7.4 


34.9 ± 3.8 


3.4 ± 0.3 


2.9 ± 1.2 


18.1 ±2.0 


4.9 ± 


1.6 


25.1 ±7.7 


2.9 


± 


0.6 


1.9±0.2 


16.0± 0.6 


83.80 


-4.25 


59 


7.3 ± 1.4 


14.4 ± 2.5 


11.2 ± 1.7 


1.1 ±0.2 


1.3 ±0.9 


13.9 ± 1.5 


3.2 ± 


0.6 


1.9 ±2.3 


2.5 


± 


0.6 


1.2 ±0.3 


16.2 ±0.5 


89.93 


2.01 


111 


2.1 ±0.3 


5.5 ±0.5 


5.7 ±0.3 


0.8 ± 0.0 


3.0 ±0.1 


7.4 ± 0.2 


1.7 ± 


0.1 


4.8 ± 0.3 


2.2 


± 


0.5 


1.5 ±0.1 


15.0 ±0.5 


93.00 


-0.07 


111 


5.0 ±0.8 


11.8 ± 1.3 


10.0 ± 0.9 


1.3 ±0.1 


3.8 ± 0.4 


8.5 ± 0.7 


1.8 ± 


0.3 


1.8 ± 1.1 


2.4 


± 


0.5 


1.8 ±0.2 


15.6 ±0.5 


99.14 


-4.25 


91 


7.6 ± 1.1 


19.5 ±2.6 


17.4 ±2.1 


0.7 ± 0.7 


0.9 ± 2.0 


20.4 ± 2.6 


<4.8 




<12.5 


4.0 


± 


0.5 




11.9 ±0.3 


99.14 


-0.07 


110 


3.3±0.9 


7.2 ± 1.3 


5.9 ± 1.0 


0.4 ± 0.2 


1.1 ±0.4 


3.1 ±0.8 


0.9 ± 


0.4 


4.8 ± 1.8 


3.1 


+ 


0.7 


1.4 ±0.5 


13.8 ±0.4 


99.14 


12.46 


79 


1.6 ±0.4 


5.7 ±0.8 


6.0 ± 0.7 


0.9 ±0.1 


2.8 ± 0.4 


7.8 ± 0.9 


1.3 ± 


0.6 


9.6 ± 3.6 


2.6 


± 


0.6 


1.6 ±0.3 


13.2 ±0.4 


102.21 


10.37 


99 


2.0 ± 0.4 


6.2 ±0.9 


6.3 ±0.7 


0.6 ±0.1 


2.9 ± 0.3 


7.0 ± 0.5 


1.0 ± 


0.3 


1.9 ± 1.0 


2.9 


± 


0.6 


1.5 ±0.2 


13.0 ±0.4 


102.21 


12.46 


68 


1.8 ±0.3 


5.8 ±0.7 


5.9 ±0.5 


0.6 ±0.1 


3.1 ±0.3 


6.4 ± 0.4 


0.3 ± 


0.3 


<3.8 


3.0 


± 


0.5 


1.5 ±0.2 


12.8 ±0.4 


108.34 


-4.25 


90 


2.6 ±0.3 


7.9 ± 0.6 


8.1 ±0.5 


1.4 ±0.1 


3.8 ±0.2 


8.9 ± 0.3 


1.7 ± 


0.1 


9.0 ± 0.8 


2.2 


± 


0.5 


1.9 ±0.1 


14.6 ± 0.5 


108.34 


2.01 


111 


15.7 ± 1.1 


24.2 ± 1.4 


14.3 ± 0.8 


1.1 ±0.1 


4.4 ± 0.2 


9.4 ± 0.5 


2.0 ± 


0.2 


1.8 ±0.4 


2.8 


± 


0.4 


1.7 ±0.1 


17.7 ±0.6 


111.41 


6.19 


83 


3.4 ± 1.6 


8.5 ± 2.7 


7.5 ± 1.5 


0.8 ± 0.2 


3.9 ± 0.4 


8.6 ± 0.7 


1.8 ± 


0.3 


4.8 ± 0.8 


2.7 


+ 


0.9 


1.7 ±0.2 


14.3 ± 0.4 


114.48 


-0.07 


110 


18.9 ± 1.7 


31.0 ±2.5 


19.5 ± 1.4 


0.8 ± 0.2 


3.9 ± 0.4 


8.0 ± 0.7 


1.9 ± 


0.2 


8.0 ± 1.1 


3.4 


+ 


0.4 


1.6 ± 0.2 


14.6 ± 0.4 
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Long. Lat. pix. bioo 


bl40 


b240 


b550 


bs50 


bl382 


b2098 


b3191 




osubmm 

Pco 


T co 


±Abioo 


±Abi4 


±Ab 24 o 


±Ab 55 o 


±Ab 85 o 


±Abi382 


+Ab2098 


±Ab3i9i 




■ a osubmm 


±ATj° 










(xlO- 1 ) 


(xl(T 2 ) 


(xlO- 2 ) 


(xlO" 3 ) 





117.55 


-0.07 


110 


1.5 


± 


1.3 


6.1 +2.6 


6.2 ± 1.7 


1.0 


± 


0.3 


2.3 


± 


0.6 


7.3 ± 0.9 


2.1 


± 


0.2 


2.4 ± 1.3 


3.2 ± 1.7 


1.4 ±0.3 


12.5 ±0.4 


123.68 


2.01 


111 


1.5 


± 


0.2 


5.2 + 0.4 


6.1 ±0.3 


0.8 


± 


0.1 


2.9 


± 


0.1 


7.3 ±0.3 


1.5 


± 


0.1 


5.1 ±0.6 


2.7 ± 0.5 


1.6 ±0.1 


13.0 ±0.4 


129.82 


2.01 


110 


3.5 


+ 


0.5 


10.8 + 0.9 


10.3 ± 0.6 


1.3 


± 


0.1 


5.7 


± 


0.3 


13.0 ±0.6 


2.4 


± 


0.3 


<4.3 


2.7 ± 0.5 


1.6 ±0.1 


13.6 ±0.4 


135.95 


-4.25 


49 


1.2 


± 


0.5 


5.1 ± 1.1 


5.4 ± 1.0 


1.3 


± 


0.2 


3.1 


± 


0.5 


10.3 ± 1.1 


3.5 


± 


0.7 


19.7 ±4.0 


2.7 ±0.9 


1.4 ±0.2 


12.9 ±0.4 


135.95 


-2.16 


106 


5.1 


± 


0.7 


11.0 ± 1.3 


8.6 ± 0.9 


0.7 


± 


0.1 


3.9 


± 


0.5 


8.0 ± 0.8 


1.3 


± 


0.4 


5.7 ± 1.4 


2.9 ±0.5 


1.6 ±0.3 


14.6± 0.4 


139.02 


-4.25 


62 


1.8 


± 


0.3 


6.6 + 0.7 


7.2 ± 0.6 


1.2 


± 


0.1 


4.0 


± 


0.3 


11.0 ±0.7 


3.2 


± 


0.3 


14.8 ± 1.2 


2.5 ±0.5 


1.3 ±0.2 


13.5 ±0.4 


139.02 


-0.07 


111 


8.6 


+ 


3.0 


21.0 + 4.1 


16.6 ±2.2 


2.6 


± 


0.3 


6.7 


± 


0.8 


14.5 ± 1.2 


3.8 


+ 


0.4 


7.3 ± 1.2 


2.4 ± 0.7 


1.8±0.2 


16.0 ± 0.6 


139.02 


2.01 


110 


6.8 


+ 


2.8 


17.3 ± 3.8 


14.8 ±2.0 


2.4 


± 


0.2 


5.7 


± 


0.6 


15.8 ±0.9 


5.0 


+ 


0.3 


14.8 ± 1.3 


2.3 ±0.8 


1.5 ±0.1 


15.8±0.5 


145.16 


4.10 


111 


2.6 


+ 


0.6 


6.7 ± 1.0 


7.0 ± 0.7 


1.1 


± 


0.2 


2.6 


± 


0.3 


4.9 ± 1.1 


1.7 


± 


0.5 


4.8 ±3.0 


2.0 ±0.6 


1.9 ±0.3 


15.5 ±0.5 


145.16 


10.37 


116 


2.1 


± 


0.1 


6.9 + 0.4 


7.3 ±0.3 


1.0 


± 


0.1 


3.1 


± 


0.2 


8.5 ± 0.5 


2.7 


± 


0.3 


8.2 ± 1.2 


2.6 ±0.4 


1.3 ±0.1 


13.4 ±0.4 


151.29 


4.10 


110 


1.6 


± 


0.2 


5.5 + 0.4 


5.7 ±0.3 


0.9 


± 


0.1 


2.6 


± 


0.2 


7.9 ± 0.4 


1.1 


+ 


0.2 


4.2 ± 0.8 


2.6± 0.5 


1.6 ±0.1 


13.5 ±0.4 


154.36 


4.10 


111 


1.8 


± 


0.2 


5.2 + 0.5 


5.6 ± 0.4 


0.8 


± 


0.1 


2.7 


± 


0.1 


7.1 ±0.4 


1.7 


+ 


0.2 


2.6 ± 0.5 


2.3 ±0.5 


1.5± 0.1 


14.4 ± 0.4 


157.43 


-0.07 


110 


1.3 


± 


0.2 


4.5 ± 0.6 


5.3 ±0.6 


1.3 


± 


0.1 


3.5 


± 


0.6 


6.5 ± 1.2 


1.6 


+ 


0.5 


5.4 ± 1.8 


2.0 ± 0.6 


2.0 ± 0.2 


14.5 ±0.5 


160.50 


-20.96 


72 


1.4 


+ 


0.2 


4.1 ±0.5 


4.5 ± 0.3 


0.5 


± 


0.0 


1.9 


± 


0.2 


5.4 ±0.3 


1.2 


± 


0.2 


4.9 ± 0.4 


2.6 ±0.5 


1.4 ±0.2 


13.5 ±0.4 


160.50 


-18.87 


91 


2.2 


± 


0.5 


5.5 ± 1.0 


5.4 ± 0.6 


0.4 


± 


0.1 


2.2 


± 


0.2 


6.0 ± 0.5 


1.1 


± 


0.1 


4.1 ±0.4 


2.9 ±0.6 


1.4 ±0.2 


13.5 ±0.4 


160.50 


-0.07 


117 


2.7 


+ 


0.3 


8.4 ±0.8 


9.3 ±0.8 


1.7 


+ 


0.1 


3.1 


± 


0.6 


6.8 ± 1.2 


2.6 


± 


0.5 


1.8 ± 1.5 


2.1 ±0.5 


1.9 ±0.2 


14.5 ±0.5 


169.70 


-6.34 


110 


0.6 


± 


0.1 


2.7 ± 0.2 


3.6 ±0.2 


0.3 


± 


0.0 


1.8 


± 


0.1 


5.4 ± 0.2 


1.2 


± 


0.1 


6.7 ± 0.8 


6.5 ± 1.5 


1.6 ±0.2 


8.4 ± 0.4 


172.77 


-18.87 


111 


0.9 


± 


0.1 


3.3 ±0.3 


4.5 ± 0.3 


1.2 


± 


0.1 


2.9 


± 


0.2 


7.8 ± 0.2 


2.5 


± 


0.1 


4.2 ± 0.6 


1.6 ±0.5 


1.5 ±0.1 


14.7 ±0.5 


172.77 


-14.69 


110 


1.5 


± 


0.2 


6.5 ± 0.6 


9.4 ± 0.8 


1.2 


± 


0.1 


4.8 


± 


0.5 


9.4 ± 0.6 


2.3 


± 


0.2 


8.0 ± 0.8 


3.1 ±0.5 


1.8 ±0.2 


11.6 ±0.3 


172.77 


-8.43 


91 


1.1 


± 


0.1 


4.3 ±0.3 


4.9 ± 0.3 


0.8 


± 


0.1 


2.5 


± 


0.3 


7.6 ± 0.3 


1.4 


± 


0.2 


3.9 ± 1.1 


2.6 ±0.5 


1.5±0.2 


13.1 ±0.4 


172.77 


-6.34 


81 


1.2 


+ 


0.2 


4.0 ± 0.7 


5.1 ±0.6 


0.7 


+ 


0.1 


2.1 


± 


0.4 


6.0 ± 0.7 


1.8 


+ 


0.5 


6.8 ± 2.5 


2.5 ±0.6 


1.4± 0.2 


13.0 ±0.3 


175.84 


4.10 


109 


8.1 


± 


0.5 


16.0 ±0.7 


12.3 ±0.5 


0.3 


± 


0.1 


2.5 


± 


0.2 


6.6 ± 0.4 


1.9 


± 


0.2 


7.5 ± 0.9 


4.3 ±0.4 


1.1 ±0.2 


12.0 ±0.3 


178.91 


-23.05 


51 


3.9 


± 


0.8 


9.7 ±2.7 


12.4 ± 2.2 


1.3 


± 


0.5 


5.3 


± 


0.8 


12.6 ± 1.9 


4.4 


± 


1.4 


12.1 ±7.6 


2.2 ± 0.6 


1.3 ±0.4 


14.3 ± 0.4 


178.91 


2.01 


110 


3.1 


± 


0.3 


7.6 ± 0.7 


7.6 ± 0.6 


0.6 


± 


0.3 


4.3 


± 


0.5 


9.3 ± 0.9 


2.2 


+ 


0.3 


4.8 ± 0.8 


2.7 ± 0.5 


1.6± 0.2 


13.7 ± 0.4 


181.97 


-12.61 


66 


0.3 


+ 


0.2 


1.3 ±0.7 


2.0 ± 0.5 


0.4 


± 


0.1 


1.9 


± 


0.1 


4.4 ± 0.5 


1.4 


± 


0.3 


4.9 ± 1.7 


2.3 ±2.1 


1.3 ±0.2 


12.5 ±0.3 


181.97 


-8.43 


48 


1.3 


+ 


0.8 


5.1 ±2.3 


8.6 ± 2.0 


1.2 


+ 


0.3 


4.9 


± 


0.7 


9.2 ± 2.0 


2.1 


± 


0.8 


4.7 ± 3.9 




2.6 ± 1.7 


6.8 ± 1.7 


181.97 


-2.16 


111 


2.1 


± 


0.3 


6.7 ±0.8 


8.3 ±0.7 


0.5 


± 


0.2 


4.2 


± 


0.3 


9.7 ± 0.6 


2.6 


± 


0.3 


11.9 ± 1.3 


3.7 ±0.5 


1.5 ±0.2 


11.5+ 0.3 


181.97 


2.01 


111 


3.8 


± 


0.4 


8.9 ±0.8 


7.7 ± 0.6 


0.6 


± 


0.2 


4.2 


± 


0.3 


7.6 ± 0.6 


1.4 


± 


0.4 


5.6 ± 1.5 


2.9 ±0.5 


1.9 ±0.3 


14.0 ± 0.4 


185.04 


-10.52 


48 


0.2 


± 


0.2 


0.6 ± 0.8 


2.0 ± 0.6 


< 0.3 




2.0 


± 


0.2 


3.6 ± 0.6 


1.1 


± 


0.3 


3.5 ± 1.8 




2.0 ±0.3 


9.8 ± 0.2 


Median 


values 
















































2.2 ±0.4 


6.7 ± 0.8 


7.5 ± 0.7 


0.9 ±0.1 


3.1 ±0.4 


8.0± 0.7 


1.8 ±0.3 


4.9 ± 1.3 


2.6 ± 0.5 


1.6 ±0.2 


13.7 ±0.4 



Case 3 



68.46 


4.10 


34 


35.1 ± 3.4 


60.2 ± 6.7 


34.9 ± 4.4 


3.4 ± 0.5 


3.4 ± 4.2 


28.5 ± 2.8 


6.3 ± 1.8 


7.7 ± 6.4 


2.7 ±0.5 


1.5 ±0.3 


17.5 ±0.6 


77.66 


-0.07 


61 


27.2 ±2.9 


40.3 ± 3.2 


23.8± 1.3 


2.5 ±0.1 


8.1 ±0.7 


17.8 ±0.8 


5.4 ± 0.4 


6.5 ± 0.5 


2.5 ± 0.4 


1.5 ±0.1 


19.6 ±0.8 


77.66 


2.01 


107 


26.3 ±2.9 


39.4 ± 3.2 


23.8 ± 1.3 


2.6 ±0.1 


8.7 ± 0.6 


18.3 ±0.7 


5.6 ± 0.4 


7.4 ± 0.6 


2.5 ± 0.4 


1.5 ±0.1 


19.5 ±0.8 


83.80 


-0.07 


105 


42.2 ± 3.0 


62.8 ±4.1 


34.5 ± 1.9 


2.7 ±0.1 


10.7 ± 0.5 


20.9 ± 0.9 


6.7 ± 0.3 


10.0 ±0.8 


2.9 ± 0.4 


1.4 ±0.1 


18.0 ±0.7 


93.00 


-12.61 


17 


1.9 ± 1.3 


2.6±3.6 


5.8 ±2.0 


0.7 ± 0.2 


2.9 ± 1.1 


6.9 ± 1.6 


2.1 ± 1.2 


10.6 ±5.7 


0.9 ± 1.9 


1.3 ±0.5 


20.7 ± 0.8 


96.07 


-6.34 


43 


11.0 ±3.1 


20.7 ± 4.0 


16.7 ± 1.8 


1.7 ±0.2 


5.0 ± 0.3 


10.7 ± 0.7 


1.3 ±0.5 


<5.9 


2.4 ± 0.6 


1.8 ±0.2 


16.5 ± 0.5 


96.07 


-4.25 


100 


9.5 ± 1.4 


19.2 ±2.0 


15.9 ± 1.0 


2.2 ±0.1 


5.4 ± 0.4 


11.2 ± 1.0 


2.1 ±0.4 


0.4 ± 1.5 


2.2 ±0.5 


2.0 ±0.2 


17.1 ±0.6 


102.21 


-0.07 


111 


5.2 ±2.6 


10.8 ±3.9 


8.5 ± 2.2 


1.4 ±0.2 


3.6 ±0.8 


8.6 ± 1.5 


2.4 ± 0.6 


3.9 ±2.7 


2.2 ± 1.1 


1.6 ±0.3 


17.7 ±0.7 


105.27 


-0.07 


111 


4.5 ±2.5 


11.8 ± 3.6 


10.2 ± 2.0 


2.6 ± 0.3 


5.3±0.6 


11.8 ± 1.1 


2.6 ± 0.4 


5.1 ± 1.6 


2.0 ± 1.1 


2.0 ±0.2 


17.1 ±0.7 



continued on next page 



continued from previous page 



Long. 


Lat. 


pix. 


bioo 


t>140 


b240 


b550 


b850 


bl382 


b2098 


t>3191 


^co 


Pco 


1 d 








+Abioo 


±Abi4 


+Ab 24 o 


+Ab 55 o 


±Ab 85 o 


±Abi382 


+Ab2098 


+Ab 3 i9i 


±A0™ 


■ a osubmm 


±AT^° 
















(xlO- 1 ) 


(xl0~ 2 ) 


(xl0~ 2 ) 


(xlO- 3 ) 




126.75 


-2.16 


108 


4.4 ± 0.9 


7.7+ 1.3 


6.2± 0.9 


0.5 + 0.1 


2.9 ± 0.4 


6.5 ± 0.6 


0.7 ± 0.3 


<1.6 


2.5±0.6 


1.5 ±0.3 


16.3 ±0.5 


142.09 


-2.16 


111 


3.9 ± 1.2 


8.2+ 2.2 


5.7 ± 1.5 


1.4 + 0.2 


1.6 + 0.6 


3.8 ± 1.1 


0.3+0.5 


2.1 ± 1.8 


2.0 ± 0.9 


2.8 ±0.5 


19.5 ±0.9 


142.09 


12.46 


55 


2.1 +0.2 


6.8 + 0.7 


6.7 ± 0.6 


0.3 + 0.1 


2.0 + 0.7 


3.3± 1.2 


1.8+0.4 


8.8 + 1.9 


4.4 ± 0.5 


1.0 ±0.4 


11.0 ±0.2 


145.16 


-2.16 


110 


3.0 + 0.5 


6.8 + 1.0 


6.3 + 0.9 


1.5 + 0.1 


3.7 ± 0.5 


9.2 ± 1.1 


1.3 + 0.4 


1.1 + 1.1 


1.7 ±0.6 


1.9 ±0.2 


18.4 ±0.7 


148.23 


-2.16 


111 


6.0 + 0.6 


11.0+ 1.0 


9.5 + 0.8 


1.4 + 0.2 


5.2 + 0.5 


11.5 + 0.9 


2.9 ± 0.3 


8.1 + 1.3 


2.0 ± 0.5 


1.5 ±0.2 


18.4 ±0.7 


154.36 


-8.43 


53 


0.9 + 0.1 


3.2 + 0.2 


4.1 +0.3 


0.5 + 0.0 


1.9+0.2 


7.2 + 0.3 


1.4 + 0.2 


4.4 + 0.7 


2.7 ±0.5 


1.2 ±0.1 


12.6 ±0.3 


154.36 


-4.25 


98 


1.1 +0.2 


4.0 + 0.6 


5.3 + 0.6 


1.3 + 0.2 


2.7 ± 0.4 


7.1 ±0.6 


1.5 + 0.4 


3.9 ± 1.8 


1.9 ±0.6 


1.9 ±0.2 


14.1 ±0.4 


154.36 


14.54 


20 


2.9 + 0.4 


11.6+ 2.0 


9.9+ 1.6 


1.9 + 0.4 


2.7 + 1.0 


11.4+ 1.9 


2.3 + 1.8 


4.3 + 6.2 


3.1 ±0.7 


1.9 ±0.4 


12.8 ±0.5 


157.43 


2.01 


111 


2.3+0.2 


5.4 + 0.3 


5.5 + 0.3 


1.1 +0.0 


3.5 + 0.2 


6.6 + 0.4 


1.2 + 0.2 


3.2 + 0.5 


1.8 ±0.4 


1.9 ±0.1 


17.3 ±0.6 


157.43 


6.19 


73 


2.3+0.1 


7.0 + 0.3 


7.4 + 0.3 


1.3 + 0.0 


4.8 + 0.2 


11.0 + 0.3 


3.0 + 0.1 


10.1 ±0.8 


2.1 ±0.4 


1.5 ±0.1 


14.7 ± 0.5 


166.63 


-2.16 


111 


2.7 + 0.3 


6.7+0.7 


7.2 + 0.6 


2.0 + 0.2 


3.6 + 0.5 


7.9 + 0.6 


1.2 + 0.3 


<3.5 


1.5 ±0.5 


2.3 ±0.2 


18.0 ±0.7 


166.63 


2.01 


111 


7.5 + 1.2 


16.0 + 3.2 


16.6 + 2.4 


3.9 + 0.5 


9.9 + 1.6 


21.3 + 2.9 


7.4 ± 1.2 


11.9 ± 5.3 


1.5 ± 0.6 


1.6 ± 0.2 


19.0 ± 0.7 


166.63 


4.10 


90 


3.0+1.0 


4.4 ± 2.6 


8.5 + 2.0 


1.9 + 0.6 


6.6+ 1.4 


17.5 + 2.5 


6.4 ± 1.7 


17.2 ±5.8 


0.6 ± 0.9 


1.1 ±0.3 


23.2 ± 1.0 


1 69 70 


-12.61 


1 10 

1 1 V./ 


1.4 + 0.1 


5 + 03 


6 2+03 


13+00 


3 6+01 


7 + 02 


13 + 01 


3 8 + 05 


2.0 + 0.4 


2.0 + 0.1 


14.1 + 0.4 


178.91 


-6.34 


93 


3.4 + 0.3 


10.6+0.7 


12.0 + 0.6 


1.0 + 0.1 


6.6 + 0.5 


10.9 + 0.5 


1.9 + 0.3 


0.7+1.2 


3.1+0.4 


1.8 ±0.2 


12.4 ± 0.3 


178.91 


6.19 


88 


2.2 + 0.3 


5.3 + 0.7 


5.4 + 0.5 


0.9 + 0.1 


4.6 + 0.3 


7.1+0.6 


0.9 + 0.3 


0.3 ± 1.1 


2.0 ±0.5 


1.8±0.2 


16.1 ±0.5 


Median 


values 
































3.4 + 0.9 


8.2 + 2.0 


8.5 + 1.0 


1.4+ 0.1 


3.7 + 0.5 


10.7 ± 0.9 


2.1 + 0.4 


4.3 ± 1.5 


2.1 ±0.5 


1.6 ± 0.2 


17.3 ±0.6 



Table 8. Results of the 1R/CO correlations in the molecular phase for all sky regions. Cols 1-2: central position of each studied regions in degree. Col 3: number of pixels used 
for the correlations in each studied regions. Cols 4-11: results of the correlations in MJy/sr for 10 21 H/cm 2 . Cols 12-13: emissivity spectral indices in the range 100 - 240 jim 
(J3 FlR ) and in the range 550 - 2100 ^m ( y g submm ). Col 14: derived dust temperature in Kelvin. Errors bars represent the 1-cr uncertainty estimates. Uncertainties equal to 0.0 are 
inferior to 0.05. 
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